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ABSTRACT 

The  classical  n-port  resistance  network  design  problem 
is  defined  and  a  common  approach  to  its  solution  based  on 
parameter  optimization  is  offered.   An  alternative  approach 
is  then  proposed  in  which  dependent  sources  are  introduced 
into  the  network  and  where  dependent  source  parameters 
assume  the  primary  roles  in  the  parameter  optimization 
problem.   Development  of  the  alternative  approach  leads 
to  specification  of  a  related  adjoint  network  and  ultimately 
to  the  generation  of  an  iterative  design  algorithm. 
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I.   INTRODUCTION 

A.   CLASSICAL  N-PORT  DESIGN 

Design  of  the  n-port  resistance  network  is  a  classical 
problem  in  circuit  theory  [1].   Necessary  and  sufficient 
conditions  are  sought  which  guarantee  realizability  of  a 
given  real,  symmetric,  nonnegative  definite  nxn  matrix  as 
an  admittance  or  an  impedance  matrix  of  an  n-port  network. 
The  network  is  to  consist  entirely  of  real,  positive-valued, 
two-terminal  resistance  elements.   Such  necessary  and 
sufficient  conditions  do  not  exist,  in  general. 

Solution  of  the  problem  using  parameter  optimization 
theory  with  the  aid  of  the  digital  computer  has  been  shown 
to  be  a  sound  practical  approach  [2].   The  parameters  to 
be  optimized  are  simply  the  individual  branch  resistances 
of  the  n-port  network.   In  this  context,  optimizing  param- 
eters implies  adjusting  the  branch  resistance  (conductance) 
values  of  the  n-port  network  to  achieve  the  desired  network 
specified  by  the  given  impedance  (admittance)  matrix.   In 
lieu  of  the  set  of  necessary  and  sufficient  conditions 
originally  sought,  application  of  parameter  optimization 
theory  and  the  subsequent  generation  of  an  iterative  design 
algorithm  result  in  convergence  of  the  algorithm  for  reali- 
zable matrices  and  no  convergence  for  unrealizable  matrices. 


B.   AN  ALTERNATIVE  APPROACH 

An  alternative  method  is  now  proposed  for  solution  of 
the  parameter  optimization  problem  where  individual  branch 
resistances  are  held  constant  and  dependent  sources  are 
introduced  into  the  network.   The  four  basic  types  of  de- 
pendent sources  are  shown  in  Fig.  1.1.   The  parameters  of 
the  dependent  sources  now  assume  the  primary  roles  in  the 
optimization  problem.   In  addition,  introduction  of  the 
dependent  sources  implies  that  more  than  one  branch  may 
now  connect  a  pair  of  nodes.   Whereas  the  original  design 
algorithm  relied  on  the  single  parameter  influence  of 
branch  resistance  variation,  the  proposed  alternative 
method  allows  multi-parameter  influence  in  the  variation 
of  each  of  the  dependent  source  parameters. 

Consistent  with  the  original  formulation,  an  exact 
solution  to  the  resistance  n-port  design  problem  safely 
may  be  assumed  to  be  impossible,  in  general.   Therefore, 
a  quadratic  error  criterion  is  employed  to  define  the 
parameter  optimization  problem.   Specifically,  the  weighted 
sum  of  the  squared  errors  between  actual  and  desired  entries 
in  the  given,  nxn  matrix  at  least  assures  convergence  to 
a  local  minimum  for  realizable  matrices.   A  typical  quadratic 
performance  criterion  would  be 
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Fig.  1.1.   Assumed  configurations  for  the  four  basic 
types  of  dependent  sources. 

(a).  Voltage-controlled  voltage  sources. 
Superscript  W  identifies  branches 
of  these  sources. 

(b).  Voltage-controlled  current  sources. 
Superscript  VI  identifies  branches 
of  these  sources. 

(c).  Current-controlled  voltage  sources. 
Superscript  IV  identifies  branches 
of  these  sources. 

(d).  Current-controlled  current  sources. 
Superscript  II  identifies  branches 
of  these  sources. 

The  C  subscript  identifies  the  controlling 
branch  and  the  D  subscript  identifies  the 
dependent  branch. 


where  p_  represnets  a  general  vector  of  parameters  to  be 


optimized  for  an  n-port  network,  and  y,   represents  the 
desired  value  of  the  k£   entry  of  a  general  nxn  admittance 
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whose  k£   entry  is  yk2« 

The  w,   are  real,  nonnegative  weighting  factors  which 
add  design  flexibility  to  the  problem.   For  example,  a 
given  entry  in  the  nxn  admittance  or  impedance  matrix  may 
be  ignored  by  choosing  the  corresponding  w,   weighting 
factor  to  be  zero.   In  addition  for  situations  in  which 
the  given  matrix  is  nonrealizable ,  appropriate  selection 
of  the  weighting  factors  is  equivalent  to  effecting  engineer- 
ing compromises  in  the  design.   Finally,  when  dealing  with 
an  iterative  design  scheme,  adjustment  of  the  weighting 
factors  during  the  course  of  iteration  may  serve  to  speed 
convergence . 


II.   FORMULATION  OF  THE  PROBLEM 

A.   N-PORT  DESIGN  AS  A  SIGNAL  APPROXIMATION  PROBLEM 

The  problem  formulation  parallels  that  of  Director  [3] 
with  modifications  for  the  introduction  of  dependent  sources 
Only  the  realization  of  a  given,  real,  nxn  admittance 
matrix  is  considered,  recognizing  the  close  similarity 
of  the  impedance  realization  problem. 

The  resistance  n-port  design  problem  with  dependent 
sources  is  treated  as  a  signal  approximation  problem.   As 
shown  in  Fig.  2.1,  the  general  network  under  consideration, 
denoted  as  71  ,  is  assumed  to  have  voltage  excitation  and 
current  response.   Since  the  problem  deals  strictly  with 
linear  time-invariant  resistance  networks  with  the  dependent 
sources,  consideration  need  be  given  only  to  n  linearly 
independent  n-vector  dc  port  voltage  excitation  signals  and 
their  n  corresponding  n-vector  dc  current  response  signals. 
The  n-port  resistance  network  with  dependent  sources  is  to 
be  characterized  by  a  given,  real,  constant,  nxn  admittance 
matrix  Y  requiring  that  the  dc  independent  voltage  source 
port  excitation  n-vector 


v  =  (vx  v2  .  .  .  vn)T 


and  the  dc  port  current  response  n-vector 


i  =  (i1  i2  .  .  .  in)T 
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Fig.  2.1.   General  n-port  network. 
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be  related  by 

i  =  Yv  .  (2-1) 

In  addition,  define  the  nxn   nonsingular  voltage  excita- 
tion matrix 


1   =  C  ^1 1  —2  I  *  *  *  l-n-1 


where 


vx  =  (vn  v21  .  .  .  vnl)T 
v2  =  (v12  v22  .  .  .  vn2)T 


T 
v   =  (v,   V~   ...  V   ) 
-n     In   2n        nn 


The  n  columns  of  V,  having  been  denoted  as  v. (j  =  l,2,...,n), 

J 

represent  a  set  of  linearly  independent  port  voltage  exci- 
tation signals.   They  define  the  corresponding  nxn  port 
current  response  matrix 


1=  [  i -,_  I  i 2 1  •  •  -I  in] 


where 


T 
-1  =  (ill  i21  *  '  *  inl) 

-2    U12  X22  *  *  *  1n2; 


T 
i   =  (i,   in   .  .  .  i   ) 

— n     In   2n        nn 
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The  more  general  form  of  equation  (2-1)  can  now  be  written 

I  =  YV  (2-2) 

Since  V  is  nonsingular  (because  its  columns  are  chosen  to 
be  linearly  independent),  its  inverse  exists  and  equation 
(2-2)  can  be  written  as 

Y  =  IV"1  (2-2a) 

where  the  choice  of  the  n  linearly  independent  dc  n-vector 
excitation  signals  is  completely  arbitrary.   For  convenience, 
therefore,  let 

V  =  U 

where  U  is  the  nxn  unit  matrix,  allowing  the  nxn  admittance 

matrix  Y  to  be  completely  specified  in  terms  of  port  current 
response  signals  by 

1=1.  (2-3) 

The  n-port  resistance  design  problem  with  dependent 

sources  can  finally  be  stated  as  follows : 

Given  a  real,  constant,  nxn  admittance 
matrix  Y,  obtain  an  n-port  network  com- 
prised exclusively  of  real,  positive- 
valued  resistances  and  of  any  combination 
of  the  four  basic  types  of  dependent 
sources  such  that  the  nxn  nonsingular 
voltage  excitation  matrix  V  =  U  yields  the 
nxn  current  response  matrix  I  =  Y  (or 
indicate  nonrealizability  when  appropriate). 
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B.   CHARACTERISTICS  OF  THE  N-PORT  NETWORK 

1.   With  Two-Terminal  Resistance  Elements 

The  resistance  n-port  network  contains  a  maximum 

of  2n  nodes  and  a  minimum  of  n+1  nodes  where  internal  nodes 

have  been  suppressed  with  the  star-mesh  transformations  [4]. 

The  minimal  node  configuration  occurs  when  there  is  a  single 

node  common  to  all  ports.   The  maximum  number  of  internal 

two-terminal  resistance  branches,  defined  to  be  n  -n,  is 

»  g   ' 

n  -n  =  (2n-l)  +  (2n-2)  +  .  .  .  +  1 
g 

=  n(2n)-n  =  n(2n-l)  . 

Addition  of  an  independent  voltage  source  at  each  port 

2 
brings  the  total  number  of  branches  to  n  =n(2n-l)+n=2n  . 

o 

In  general,  the  total  maximum  number  of  two-terminal 

branches,  n  -n,  connecting  n  nodes  is 
9   g   *  &  n 

n  -n  =  (n  )!/2!(n  -2)!  =  n  (n  -D/2. 
g       n       n        n  n 

The  minimal  branch,  minimal  node  configuration  is  found 
to  be 

n  -n  =  ((n+l)-l)  +  ((n+l)-2)  +  .  .  .  +  1 

2 
=  n(n+l)-n  =  n  -n  . 

Again,  addition  of  an  independent  voltage  source  at  each 

2        2 
port  brings  the  total  number  of  branches  to  n  -n+n  =  n  . 
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Summarizing,  the  general  n-port'  resistance  network 

contains  n  nodes,  where  (n+l<n  <2n),  and  n  total  internal 
n      '  —  n—         g 

2      2 
plus  excitation  branches,  where  (n  <n  <2n  ). 

2 .   With  Two-Terminal  Resistance  Elements  and 

Dependent  Sources 

The  significance  of  the  addition  of  dependent 
sources  in  the  problem  formulation  lies  in  the  fact  that 
more  than  one  branch  may  now  connect  a  given  pair  of  nodes 
in  the  network.   Moreover,  each  node-pair  which  is  connected 
by  a  single  two-terminal  resistance  element  is  allowed  the 
influence  of  every  other  similar  node-pair  in  the  network 
by  using  the  appropriate  dependent  sources. 

A  typical  pair  of  nodes  is  shown  in  Fig.  2.2a  which 
illustrates  the  introduction  of  dependent  sources  into  an 
n-port  resistance  network.   Each  summation  set  equal  to 
zero  represents  a  set  of  D  dependent  source  controlling 
branches,  where  D  is  equal  to  the  number  of • internal  re- 
sistance branches  in  the  network.   For  example,  ZVP   =  0 

D  ° 
represents  a  set  of  D  current-controlled  voltage  source 

controlling  branches.   Similarly,  the  remaining  summations 
each  represent  a  set  of  C  dependent  source  dependent 

branches,  where  C  is  also  equal  to  the  number  of  internal 

W 
resistance  branches  in  the  network.   For  example,  E  Vn 

C  u 

represents  a  set  of  C  voltage-controlled  voltage  source 

dependent  branches. 

■■  "     'Accounting  for  each  branch  in  the  network  is  critical 

in  the  subsequent  derivation  of  the  gradient  expression, 


14 


I   I 


II 

D 


I    I 


VI 
D 


Fig.  2.2a.   General  node-pair  configuration 
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Fig.  2.2b.   Simplified  node-pair  configuration 
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therefore  a  standard  notation  is  now  introduced  for  numbering 

the  network  branches.   First,  the  port  excitation  branches 

are  numbered  from  1  to  n.   The  resistance  branches  are  then 

numbered  from  n+1  to  n  .   Finally  the  dependent  sources  are 

g 

numbered  from  n  +1  to  n,  ,  where  the  controlling  branch  and 

g       b* 

the  dependent  branch  of  each  dependent  source  are  numbered 

sequentially,  in  that  order.   For  example,  the  controlling 

branch  of  the  first  dependent  source  considered  is  numbered 

n  +1  while  the  dependent  branch  is  numbered  n  +2. 
g  g 

All  that  remains  in  the  problem  formulation  is  to  specify 

the  limits  on  the  total  number  of  nodes  n  and  of  branches 

n 

n.  due  to  the  introduction  of  dependent  sources  into  the 

network.   Referring  to  Fig.  2.2a,  each  summation  in  series 

with  the  branch  resistance  represents  n  -n  possible  voltage 

g 

sources.   Since  there  are  four  different  summations,  the 

maximum  number  of  nodes  n  in  the  network  is  found  to  be 

n 

n   =  2n  +  4(n  -n)  =  kn   -2n. 
n  g        g 

Likewise,  each  of  the  four  different  summations  connected 
to  node-pair  AB  represents  n  -n  possible  current  sources. 
The  maximum  total  number  of  internal  plus  excitation 
branches  is  found  to  be 

n.  =  n  +  4(n  -n)  =  5n  -*ln. 
b    g      g        g 

To  simplify  the  topological  analysis  of  the  network, 
only  the  nonzero  branches  of  Fig.  2.2a  need  be  considered, 
reducing  the  maximum  number  of  nodes  and  branches  required 
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to  define  the  network.   Accordingly,  Fig.  2.2b  illustrates 
the  application  of  superposition  to  the  general  node-pair 
configuration  of  Pig.  2.2a.   The  resulting  simplified  node- 
pair  configuration  leads  to  an  equivalent  but  simplified 
form  of  the  original  network. 

Referring  to  Fig.  2.2b,  the  presence  of  one  or  more 
voltage-controlled  voltage  source  dependent  branches  or  of 
one  or  more  current-controlled  current  source  dependent 
branches  creates  an  additional  node  (denoted  B')  for  each 

of  the  existing  n  -n  resistance  branches.   The  maximum 

g 

number  of  nodes  becomes  2n+(n  -n)  =  n  +n.   Similarly,  the 

presence  of  one  or  more  current-controlled  current  source 

dependent  branches  or  of  one  or  more  voltage-controlled 

current  source  dependent  branches  creates  an  additional 

branch  which  also  connects  none-pair  AB.   Therefore,  the 

maximum  number  of  resistance  branches,  simplified  voltage 

source  dependent  branches,  and  simplified  current  source 

dependent  branches,  denoted  by  n  -n,  is 

s 

n  -n  =  n  -n+2(n  -n)  =  3(n  -n)  . 
s      g      g         g 

An  additional  excitation  branch  at  each  port  brings  the 
total  maximum  number  of  branches  for  the  topological  analy- 
sis to 

n„  =  3(n  -n)+n  =  3n  -2n  . 
s      g  g 

For  convenience  in  the  simplified  network,  excitation 
branches  are  numbered  from  1  to  n,  resistance  branches  are 
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numbered  from  n+1  to  n  ,  and  dependent  source  dependent 

branches  are  numbered  from  n  +1  to  n  . 

g       s 

Note  that  zero-valued  voltage  and  current  sources  as 
shown  in  Fig.  1.1  are  indeed  legitimate  branches  in  the 
original  unsimplified  network.   Their  inclusion  in  the  prob- 
lem formulation  will  become  more  clear  when  the  adjoint 
network  is  formally  defined  in  the  next  section. 

Summarizing,  the  general  n-port  resistance  network  with 

dependent  sources  contains  n   nodes  where  (n+l)<n  <(4n  -2n) , 

2 

and  n,  branches  where  (n  )<n,  <(5n  -4n) .   The  simplified 
b  —  b—   g  * 

form  of  the  network  contains  n  nodes  where  (n+l)<n  <(n  +n) . 

n  —  n—  g    5 

2 
and  n  branches  where  (n  )<n  < ( 3n  -2n). 
s  —  s—   g 

Complete  specification,  therefore,  of  the  given  nxn 
admittance  matrix  Y  assumes  n  separate  measurement  situa- 
tions for  the  n-port  resistance  network  with  dependent 
sources.   Each  measurement  situation  claims  a  set  of  n  dc 
port  voltage  excitation  signals  and  n  dc  port  current 
response  signals  and  an  overall  set  of  n,  dc  branch  voltage 
and  current  signals. 
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III.   DERIVATION  OF  THE  DESIGN  SCHEME 

A.   OBTAINING  THE  GRADIENT  EXPRESSION 

The  assumption  has  been  made  that  the  given  nxn  admit- 
tance matrix  may  not  be  realized  exactly.   In  addition, 
since  the  admittance  matrix  is  completely  specified  in  terms 
of  network  port  current  response  signals  by  equation  (2-3), 
Y.  =  I_,  an  approximation  criterion  may  be  written  from 
equation  (1-1)  as 


J(E)  ='  |  .?  ?  "w^ia-ikJ2  (3-D 

l-l  k=l 


where  p  now  represents  the  ^(n  -n)-vector  of  dependent 

o 

source  parameters  to  be  optimized.   As  before,  the  w,    are 

the  weighting  factors,  and  i,   is  the  desired  value  and 

i,   is  the  actual  value  of  the  k   dc  port  response  current 

due  to  the  l        linearly  independent  dc  n-vector  port  voltage 

excitation. 

A  relationship  is  sought  between  dependent  source 
parameters  and  the  given  squared  error  approximation  or 
performance  criterion.   This  knowledge  should  allow  adjust- 
ment of  the  parameters  such  that  J(p)  can  be  minimized. 
Clearly,  the  gradient  of  J(p)  with  respect  to  dependent 
source  parameters  p  is  the  desired  relation.   Director  and 
Rohrer  [2]  have  outlined  the  derivation  of  the  gradient 
expression  for  the  pure  resistance  n-port  network  upon 
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application  of  Tellegen's  theorem  [5].   The  derivation 
which  follows  is  based  on  methods  of  variational  calculus, 
and  it  parallels  that  of  Director  [31. 

Minimization  of  the  performance  measure  is  subject  to 
constraints  imposed  by  network  topology  and  by  branch 
relations.   A  topological  constraint  with  respect  to 
Kirchhoff's  current  lav;  is  written  as 


nb 

Z      q  ,  i,  =0    m  =  1,2, . . . ,n  -1 
.  ,  Mmk  k  '  5    '  n 


where  Q-=  [q  .  ]  is  any  n  -1  by  n,  fundamental  cutset  matrix, 
and  i   is  the  k   branch  current.   Kirchhoff's  voltage  law 
specifies  a  second  topological  constraint  which  is  written 
as 


nb 

I     b  .  v,  =  0    m  =  l,2,...,n,-n  +1 
.  ,   mk  k  3    3         3   b  n 


where  B  =  [b  .  ]  is  any  n,  -n  +1  by  n   fundamental  loop  matrix 
and  v,  is  the  k   branch  voltage. 

Branch  relations  for  the  network  elements  are  written  as 


v,  -v,  =  0         for  k   independent  voltage 

source,  where  v,  is  the  desired 

'        k 

value  and  v,  is  the  actual  value; 
k 

G,  v,  -i   =  0       for  k   conductance  branch; 
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C-° 


for  voltage-controlled 


.-,,    yy  voltage  sources; 


v-  -  uv--   0 


31-" 


for  voltage-controlled 
,,T    VT  current  sources; 

I   -E  V   =  0 


vf  =  0 


for  current-controlled 


TV    TV  voltage  sources; 

V:!:v-r  lt,y    =  0 
D   m  C 


i1-" 


I?-BI»  =  0 


for  current-controlled 
current  sources. 


Consistent  with  the  notation  in  Fig.  1.1,  the  C  and  D  sub- 
scripts indicate  controlling  and  dependent  branches, 
respectively,  for  the  dependent  sources. 

The  topological  constraints  and  branch  relations  are 
now  appended  to  the  performance  measure  using  Lagrange 
multipliers.   The  augmented  performance  measure  is  then 
written  as 
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,    n  p 


£=1      k=l 


n  -1  n,           n,  -n  +1  n, 

n  .    b            b   n  b 

+   E  X  „      E  q  ,  i.  „   +     E     X  .  E   b  ,,v,  „ 

.  ml   .  ,  ^mk  kJl        n     mi  ,  ,   mk  k£ 

m=l  k=l            m=l  k=l 


n  n 

+   E  xXn(vln-V1n)   +     E    X?n  (G.  v,   -i,  „  ) 

,  n   k£v  k£   k£      .   ...   k£   k  k£   k£ 
k=l  k=n+l 


n,_ 

+    .       *  [Xk£l    Jk*    +    Xk£2    ^kA-VW 

k=ng+l 


.VI      TV1         ,VI       ,TVI  ,rVIv 

+    Xk£l    Jk*    +    Xk£2    ^rk*-«nkVkA) 


,IV      ,,IV         ,IV       ,„IV  TIVv 

+    XkU   Vkl    +    \l2    (Vki-rmkIki) 


+  *k*iVk*+  ^2  Clg-3k^>3>-    (3-2) 


Although  the  network  parameters  doe  not  appear  explicitly 
in  the  original  performance  measure  (3-1),  their  influence 
is  revealed  in  terms  of  the  port  current  response 
i,  (k=l,  2,  .  .  . ,  n) .   Hence,  if  a  given  network  parameter 
is  incrementally  varied 
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p    -*■  p  +  6p 
then 


and 


\  -  \  +  61k 


vk  *  \  +  6v 


Upon  introducing  an  incremental  variation  in  all  dependent 
source  parameters,  branch  relations  for  the  network  elements 
become 

(v,+6vk)  -  v.  =  0  for  k   independent  voltage 

source; 

G,  (v,  +6v,  )  -  (i,  +61.  )  =  0    for  k   conductance  branch: 
k  k   k      k   k  ' 


I„  ■'-+  61   =  0  for  voltage-controlled 

voltage  sources 

(VD  +6VD  )  -  (y+6y)(Vc  +6VC  )  =  0 


I^1  +  Sli1   =  0  for  voltage-controlled 

current  sources; 


Vp  +6v*c   =  0  for  current-controlled 

voltage  sources; 

(vf+<).(rm+6rm)(Iiv+6Iiv)  =  0 
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V^+fiV*1  =  0  for  current-controlled 

current  sources. 

(I^+fil^1)  -  (6  +  6B)(lJI+6lJI)  =  0 

Note  that  although  branch  conductances  are  held  constant, 
variations  of  conductance  branch  voltage  and  current  appear 
due  to  the  introduction  of  the  dependent  sources. 

Having  incremented  all  dependent  source  parameters, 
equation  (3-2)  can  be  rewritten  as 


1=1  k=l 

n    .    b  b   n        t> 

+  l        X     E   q  .  (i.  +61.   „)  +     E     X \       E   b  ,  (v,  +6v.  „  ) 
,  ma   ,  n  Mmk  k£   ka        n    m£  ,  ,   mkv  k£   k£ 
m=l    *  k=l  m=l       k=l 

n  n 

.    +1  kul  k£   k£     k£2   k£   k£      k   pk   k£   k£y 


,VI  ,TVI..TVIx     ,VI  //TVI.xTVIv     f  LS  \f^^j.x^^\\ 

+    Xkil(Ik*+6Iki)    +    Xk*2((Ik*+6Ik*}    "    (SV^\)(Vk^Vki)] 


,     .IV    ,VIV,.VIV.  TV    //,rIV,x,.IVx  /  .,  WTIV,„TIVss 

+    Xk£l(V  +  6Vk*}    +    Xk£2((Vk£  +  5Vkii)    ~    (rmk+6rmk)(Ik£  +  6lk£)) 


(3-3) 
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Subtracting  (3-2)  from  (3-3)  and  neglecting  second  order 
terms  leaves  the  first  order  variation  of  the  augmented 
performance  measure  written  as 

<5J(p_)    =    J(p_+5£)  -  J(p) 

n  n 

n   -1  nL  n.  -n  +1  n, 

n  b  b      n         .  b 

+       I        x\         E      q  ,,6i.       +  X         X    „         E.    b   .  6v,  . 

,       mi      ,     .    ^mk     kz  ,        mi      .     ,      mk     ki 

m=l  k=l  m=l  k=l 

n  n 

v  S      n 

+      l      X.  n6v,       +         E      X,  n  (G,  6v.     -61.  n  ) 
,    ,      kl     ki       ,       .,    k£v   k     k£        k£ 
k=l  k=n+l 


n. 

+    £    IX.  n16I.  „  +  X,  .„(6V,  „-y,  <5V.  „  -V.  „oy.  ) 
.    . ,  L  kil     ki  k£2v  ki      k     ki      ki      ky 

k=n  +l 
g 


,VI  XTVI    ,VI  ,XTVI     xirVI  ,rVIx    s 
+  Xk^l5Ik£  +  Xk£2(6Ik£-Smk6VkrVk£6Smk) 


.IV  .,rIV  ,  ,IV  /-,rIV     XTIV  TIV.    v 

+  X,  „n6V.  .  +  X.  n(6V.  -rm,  61.  -I,  „6rm.  ) 

kil     ki  ki2        ki        k     ki     ki        k'- 


+  xlli<l +  ^"ErV^i-1**6"^  ])-      (3-^ 


After  rearranging  terms,  equation  (3-*0  can  be  written 
in  a  more  useful  form  as 
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5J(p_) 


n 

E 
1  =  1 


nT 


n   -1 


nt 


n 


n 


,       raJl  .    ,    Mmk         ,     _    kH     k£     k£        .       .,    k£ 
m=l  k=l  k=l  k=n+l 


+   .       E+1(Xk)ll  +  Xk£l  +  Xkil2-rinkXkil2+  (1-6k)XkS,2)] 
k=n  +1 


n,  -n  +1  n,  n                      n 

b      n  o      -   to  v                s      G 

E.  X     „       Z  b    .     +    [       I       X,  +         I       X.     G, 

m=1  m£    k=1  mk                        k£        k=n+1   k£    k 


6i 


kz 


n, 

v  //n  N,W  ,VI       ,,IV      ,,IV      ,,II      v       n 

+    „     E+1((1-yk)Xk£2-SmkXk,2  +  Xk£l  +  Xkil2+XkU)     ] 
k=n  +1 

g 


6v 


kl 


n. 


vu       ,,W  TTW.       ..VI  TTVI.         ,,IV  -IV.         ..II    TII,Q    x 
„     Z  +1  (Xk£2Vk£6^k+Xkil2Vk^Smk+Xk£2Ik£6rmk+Xk£2Ik£63k) 

k=n_+l 


} 


(3-5) 


Equation  (3-5)  has  thus  been  written  as  a  summation  of  three 
general  terms:   (1)  terms  multiplied  by  current  variation 
6i,  ,  (2)  terms  multiplied  by  voltage  variation  6v,  ,  and 
(3)  terms  which  include  variation  of  all  dependent  source 
parameters.   Hence,  the  desired  gradient  information  lies  in 
the  third  general  term  of  equation  3-5.   The  final  form  of 
the  gradient  expression  on  which  the  iterative  design  algo- 
rithm is  to  be  based  will  appear  following  definition  and 
treatment  of  the  adjoint  network. 
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B.   DESCRIPTION  OF  THE  ADJOINT  NETWORK 

Attention  is  now  given  to  the  description  of  the 
so-called  adjoint  network  in  order  to  gain  a  useful  inter- 
pretation of  the  Lagrange  multipliers  introduced  in  the 
augmentation  of  the  performance  measure.   In  addition  to 
the  original  network  J\ ,  consider  an  adjoint  network  /^ 
which  is  topologically  equivalent .   Tree  branch  voltages 
and  network  branch  voltages  in  the  adjoint  network  (as  in 
any  electrical  network  [6])  must  be  related  by 


[q  ,  lT[v  ]  -  [v.  ]   =   [0]  (3-6) 

mk    mk 


where  v   is  an  (n  -1  by  l)-vector  of  tree  branch  voltages 
and  v,  is  an  (n  by  l)-vector  of  network  branch  voltages. 
Likewise,  link  currents  and  network  branch  currents  in  the 
adjoint  must  be  related  by 


EV]X3  "  [V   =    co:  (3-7) 


where  i   is  an  (n,  -n  +1  by  l)-vector  of  link  currents  and 

i,  is  an  (n  by  l)-vector  of  branch  currents.   Since  the 
k         n 

only  restriction  on  the  adjoint  network  /\   is  that  it  be 
topologically  equivalent  to  the  original  network  /^,  the 
branches  of  7\^   are  free  to  be  specified  in  some  convenient 
manner. 

Consider  the  first  general  term  of  equation  (3-5). 
Define  A,  as  the  k   tree  branch  voltage  in  /^.   Now  define 

A  +■  Vi  c 

-   w,  (i,  -i,  )  as  the  k   port  dc  voltage  excitation  and  X, 
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as  the  k   conductance  branch  voltage  in  // .   Similarly 


a    **  ^     ,VI     ,VI      ,IV     ,11     ,  .  ,11 

define  -  Xkl,  -  Xfel,  -  Xk2,  ™k*k2>    ~    xk2>  and  6kXk2  as  the 

remaining  branch  voltages  in  A^.   With  these  definitions, 
the  first  general  term  of  equation  (3-5)  is  seen  to  be 
identical  in  form  to  equation  (3-6),  and  this  term  vanishes. 
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Now  consider  the  second  general  term  of  equation  (3-5). 
fine  X.  as  the  k   link  current  in  7\ .   Define  X,  as  the 
current  in  the  k   port  dc  voltage  excitation  branch  and 
X,  G,  as  the  k   conductance  branch  current  in  /(_ .   Similarly 
define  Xk2,  -  yk*k2,  -  g™kxk2>  xkl,  Xk2,  and  Xkl  as  the 
remaining  branch  currents  in  7^.   With  these  definitions, 
the  second  general  term  of  equation  (3-5)  is  seen  to  be 
identical  in  form  to  equation  (3-7),  and  this  term  likewise 
vanishes . 

Having  eliminated  the  first  two  terms  of  equation  (3-5), 
it  can  now  be  written  as 


n    n. 

E  1-1   k=n +1    "        k^2Vkl6yk   +    6k*2Vk    6smk 

ng 


,IV   TIV.  ,11  'II. a       s 

+    \l2\l6rmk   +    Xk£2Jk)l6ek    }> 


or  using  vector  notation  as 


6J(p_)  =   E  {  [vj(p_)]T[6P]  }  ,  (3-8) 

£=1 
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where  the  (n.  -n  by  1)  gradient  vector  vJ(p)  is  defined  as 


VJ(p_)   = 


XW  VW 
Ak2  Vk 


VI   VI 

k2  vk 


k2   k 


Xk2  Ik 


(3-8a) 


and  the  (n.  -n  by  1)  parameter  variation  vector  SP  is 
b   g  — 


defined  as 


6P 


6y, 


6gm1 


6rm, 


66, 


(3-8b) 


The  remaining  definitions  for  the  adjoint  network  follow 
directly  from  the  derivation  and  subsequent  manipulation  of 
equation  (3-5).   To  simplify  notation,  the  l   subscript  will 
not  be  included  in  defining  branch  voltages  and  currents 
in  /(_ ,  although  it  actually  exists  for  each  of  the  n  differ- 
ent n-vector  excitation  situations  in  the  adjoint.   Now, 
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define  for 


T.  =  -w,  (i  -i.  )         k   port  dc  voltage  excitation 


k   kv  k  k 


branch  voltage; 


$,  =X,  k   port  dc  voltage  excitation 

branch  current ; 

C       C  t  hi 

¥,=X,  k   conductance  branch  voltage; 

C  C  t  h 

$  =G.  X,  k   conductance  branch  current, 
k  k  k 


For  current-controlled  current  sources,  define 

W 

4*   =0  controlling  branch  voltages; 

$   =  X„2  controlling  branch  currents; 

$P  =~yn^n  dependent  branch  currents; 

W   W 

TC  =~*C1  dependent  branch  voltages. 


For  voltage-controlled  current  sources,  define 

VI 
$_.  =0  controlling  branch  currents; 

VI   VI 
^D  ~~^C2  controlling  branch  voltages; 

VI     VI 
$c  =gmDyD  dependent  branch  currents; 

VI   VI 
^C  =  ~*n  dependent  branch  voltages. 


For  current-controlled  voltage  sources,  define 

IV 
H*  =0  controlling  branch  voltages 
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$   =App  controlling  branch  currents; 

H*  =rnL$n  dependent  branch  voltages; 

$  =Xp.|  dependent  branch  currents. 


For  voltage-controlled  voltage  sources,  define 

$n  =0  controlling  branch  currents; 

II   II 

y     =-Xrp  controlling  branch  voltages; 

II     II 

Yp  =-Br)H'n  .  dependent  branch  voltages; 

$r  =-Xr-,  dependent  branch  currents. 


In  summary,  the  adjoint  network  f\^   is  related  to  th< 

original  network  f\   in  the  following  manner. 

(i)  Yl^  and  /(.are  topologically  identical. 
(ii)   All  conductance  brances  of  value  G  in 
l\   are  associated  with  conductance 
branches  of  value  G  in  /\ . 
(iii)   Voltage-controlled  voltage  sources  with 
parameter  p  in  /[  are  associated  with 
current-controlled  current  sources  with 
parameter  -y  in  i\  ,  and  the  roles  of 
controlling  and  dependent  branches  in 
f\  are  reversed  in  /[,  . 
'  (iv)   Current-controlled  current  sources  with 
parameter  $  in  /^  are  associated  with 
voltage-controlled  voltage  sources  with 
parameter  -6  in  f[,  and  the  roles  of 
controlling  and  dependent  branches  in 
^  are  reversed  in  f[_ . 
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(v)   Voltage-controlled  current  sources  with 

parameter  g^  and  current-controlled 

voltage  sources  with  parameter  r   in  7) 

are   associated  with  voltage-controlled 

current  sources  with  parameter  g  and 

°m 

current-controlled  voltage  sources  with 
parameter  r  ,  respectively,  in  r\   with 
the  roles  of  controlling  and  dependent 
branches  in  /]_    being  reversed  in  /£.  . 

In  conclusion,  utilizing  previous  definitions,  equation 

(3-8a),  which  is  the  expression  for  the  gradient  of  the 

performance  measure  with  respect  to  dependent  source 

parameters,  can  be  written  in  final  form  as 


vj(p_) 


W  „W 

'd  vc 


,VI  vVI 
*D    C 


IV  Tiv 

'd   c 


v11  I11 


(3-9) 


In  order  to  minimize  the  performance  measure,  the  dependent 
source  parameters  are  to  be  adjusted  in  the  negative 
gradient  direction,  -  vj(p).   Physically,  the  gradient 
information  is  seen  to  be  products  of  controlling  branch 
voltages  or  currents  in  /(__  and  corresponding  dependent 
branch  voltages  or  currents  in  t\_  . 
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C.   A  NOTATIONAL  EXAMPLE 

An  example  will  now  be  given  to  illustrate  notation, 
to  illustrate  the  correspondence  between  the  original  and 
adjoint  networks,  and  to  identify  components  of  the  gradient 
expression. 

Figure  3.1  shows  a  simple  2-port  resistance  network,/?, 
with  3  different  types  of  dependent  sources.   Since  the 
example  is  intended  to  illustrate  notation,  numerical  values 
for  the  resistances  and  for  dependent  source  parameters  are 
not  specified.   Branches  are  numbered  in  accordance  with  the 
convention  previously  established  in  problem  formulation. 

Figure  3.2  shows  the  related  adjoint  network,  /£  , 
generated  with  the  aid  of  relationships  defined  in  the 
previous  section.   Note  that  the  original  and  adjoint 
networks  are  topologically  equivalent. 

The  components  of  the  desired  gradient  information  can 
now  be  identified  in  terms  of  voltages  and  currents  in  /£, 
and  in  /(_ .   The  gradient  expression,  equation  (3-9) >  for 
this  example  is 


vJ(p)  = 


*6  V6 


IV  TIV 
'8   J8 


,11  tII 
10   10 


(3-10) 
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y6V6 


Wv^ 


Fig.  3.1.   Example  2-port  network  (original). 
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10 


IV 

^  v=o 


Fig.  3.2.   Example  2-port  network  (adjoint). 


3^ 


and  the   parameter  variation   is   given  by 


6P    = 


6  y. 


<5rm< 
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The  first  term  of  equation  (3-10)  is  seen  to  be  the 
negative  of  the  product  of  the  voltage-controlled  voltage 
source  controlling  branch  voltage  in  /^  and  the  current- 
controlled  current  source  dependent  branch  current  in  /{^. 
Similarly,  the  second  term  is  seen  to  be  the  negative  of 
the  product  of  the  current-controlled  voltage  source  con- 
trolling branch  current  in  /(_  and  the  current-controlled 
voltage  source  dependent  branch  current  in /(_ .   Finally, 
the  third  term  of  equation  (3-9)  is  the  product  of  the 
current-controlled  current  source  controlling  branch 
current  in  /£  and  the  voltage-controlled  voltage  source 
dependent  branch  voltage  in  /£ . 

D.   A  PROPOSED  ITERATIVE  DESIGN  ALGORITHM 

An  iterative  algorithm  for  the  design  of  resistance 

n-port  networks  with  dependent  sources  is  now  proposed. 

(i)   Initially,  let  all  dependent  source 
parameters  be  equal  to  zero, 
(ii)   Apply  the  £th  (£=1,  2,  .  .  . ,  n)  of 
the  n  linearly  independent  n-vector 
dc  voltage  excitations  to  /^ . 
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(iii)   Measure  all  dependent  source 
controlling  branch  voltages 
and  currents  in  /(, . 

(iv)   Calculate  the  current  response  errors 
at  the  n  ports. 
(v)   Apply  these  n  port  current  response 
errors  as  voltage  excitation  at  the 
corresponding  ports  of  7]^. 

(vi)   Measure  all  dependent  source  control- 
ling branch  voltages  and  currents  in 
/^  and  multiply  them  with  correspond- 
ing dependent  source  controlling 
branch  voltages  and  currents  in  /[ . 
This  forms  the  z        set  of  the  maximum 
l\(n   -n)  components  of  the  gradient 

o 

expression, 
(vii)   Let  1  =  1+1   and  return  to  step  (ii) 

being  certain  upon  reaching  step  (vi) 
to  add  newly  calculated  components  of 
the  gradient  expression  to  the  corre- 
sponding previously  calculated  components. 
Execute  this  step  for  (1=1 3    2,  .  .  . ,  n) 
thereby  obtaining  the  complete  gradient. 
(viii)   Adjust  the  dependent  source  parameters 
in  the  negative  gradient  direction, 
(ix)   Repeat  the  entire  procedure  from 
step  (ii)  until  the  performance 
measure  has  been  suitably  minimized 
or  until  some  other  termination 
criterion  has  been  satisfied. 

A  number  of  comments  are  in  order  with  regard  to  this 

proposed  algorithm.   Step  (ii),  for  a  given  I ,  normally 

implies  inserting  a  1  volt  source  across  the  k£    port 

where  k=£  and  inserting  0  volt  sources  across  the  remaining 
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n-1  ports.   Steps  (iii)  and  (vi)  assume  the  availability  of 
a  general  circuit  analysis  program  such  as  the  IBM  ECAP 
program  [7]  for  measuring  voltages  and  currents  in  both  the 
original  and  adjoint  networks. 

A  final  comment  deals  with  step  (viii).   The  dependent 
source  parameters  are  to  be  adjusted  in  the  direction  of 
steepest  descent  for  each  iteration  of  the  design  algorithm; 
that  is 

Ri+1  =  g1  _  a1  [vJCp1)]    1=1,  2,  .  .  . 

where  superscript  i+1  or  i  indicates  the  iteration  number 
and  a  is  a  real,  nonnegative  step  size  determining  constant 
which  may  vary  from  iteration  to  iteration.   The  proper 
choice  of  a  is  critical  using  the  steepest  descent  algo- 
rithm outlined  here  and  a  good  discussion  regarding  the 
selection  of  a  and  the  use  of  the  steepest  descent  technique 
can  be  found  in  [8]. 
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IV.   CONCLUSION 

The  classical  n-port  resistance  design  problem  has  been 
defined.   A  practical  approach  to  its  solution  has  been  to 
cast  the  problem  as  a  signal  approximation  problem  and  then 
to  apply  parameter  optimization  theory  to  minimize  an 
established  approximation  criterion. 

The  alternative  approach  offered  involves  the  introduc- 
tion of  dependent  sources  into  the  network  holding  branch 
resistance  values  constant.   The  dependent  source  parameters 
then  assumed  the  primary  roles  in  the  parameter  optimization 
problem.   Derivation  of  the  gradient  expression  led  to  the 
generation  of  an  iterative  design  algorithm  which  could  be 
programmed  on  the  digital  computer. 
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